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Abstract 

We say that the polynomial sequence (Q^^^) is a semiclassical Sobolev polynomial 
sequence when it is orthogonal with respect to the inner product 

(Pi r)g = {u,pr) + X (u, ^r) , 

where u is a semiclassical linear functional, & is the differential, the difference or 
the g-difference operator, and A is a positive constant. 

In this paper we get algebraic and differential/difference properties for such poly- 
nomials as well as algebraic relations between them and the polynomial sequence 
orthogonal with respect to the semiclassical functional u. 

The main goal of this article is to give a general approach to the study of the 
polynomials orthogonal with respect to the above nonstandard inner product re- 
gardless of the type of operator & considered. Finally, we illustrate our results by 
applying them to some known families of Sobolev orthogonal polynomials as well 
as to some new ones introduced in this paper for the first time. 
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1 Introduction 



In the pioneering work [9], D.C Lewis introduced the inner products involving derivatives 
in order to obtain the least squares approximation of a function and its derivatives. Since 

* Corresponding author: e-mail: balcazar@ual.es 



Preorint submitted to Journal of ADDroximation Theorv 



Seotember 5. 2011 



that article, a large amount of papers have appeared on orthogonal polynomials with 
respect to inner products involving derivatives (and later differences or g-differences) , the 
so-called Sobolev (or A-Sobolev or g-Sobolev) orthogonal polynomials. 

We can see the evolution of this theory through different surveys from 1993 to 2006 
(ordered by year: [10], [18], [15], [16], and [11]). In the most recent surveys some directions 
for the future are proposed. 

The most studied cases correspond to the orthogonal polynomials with respect to an inner 
product of the form 

{p-,r)s = {viQ.pr) + \{\ii,^^p^^r), p,reP, (1) 

where Ug and Ui are definite positive linear functionals, and ^ is the differential operator. 
Later, in [1] the discrete case was considered where ^ = A was taken as the forward finite 
difference operator, and in [2] the g-difference operator = Qiq was considered. 

As we can observe in the surveys mentioned above, and in the references therein, all the 
results have been obtained by considering a fixed operator, that is, the differential opera- 
tor in the continuous case, the forward operator in the discrete case and the g-difference 
operator in the g-Hahn case. Notice that a great part of the results have been obtained 
from knowing an explicit representation for the linear functionals uq and ui, for example, 
when (uq, ui) is a coherent (or symmetrically coherent) pair of functionals. The coherence 
is a case of special relevance from the numerical point of view (see [6]), although according 
to the classification of the coherent pairs (see [1], [2], and [19] for each case), the coherent 
functionals are close the classical ones; the complete classification of the coherent pairs in 
the continuous, discrete and g-Hahn cases can be found in [19], [1] and [2], respectively. 
Furthermore, a lot of papers consider the case u := uq = ui when this functional is ex- 
plicitly known, for example, if u is a very classical functional (Jacobi, Laguerre, Meixner, 
etc.) or if u is semiclassical (Freud-type). For the corresponding Sobolev orthogonal poly- 
nomials algebraic, differential/ difference and asymptotic properties are obtained as well 
as a relation with the standard polynomials orthogonal with respect to u. 

In this paper we do not work in the framework of coherence (either A-coherence or 
g-coherence) in the sense of coherence worked with in [19] ([1] or [2]), although some 
families of coherent pairs of this type are included in our approach. In fact, we consider 
the case u := uq = ui in (1) where u is a semiclassical linear functional and, therefore, 
our approach is not related to the coherence case. 

We think that the theory of Sobolev orthogonal polynomials is a big puzzle in which many 
pieces have been connected during the last few decades but it lacks a general approach. 
Taking this into account, our main objective is to give a unified approach for the theory 
of Sobolev orthogonal polynomials regardless of the operator considered. More concretely, 
we consider a semiclassical linear functional u and the nonstandard inner product 

(P, '^)5 = (u,J5r) + A(u, ^p^r), rGP, (2) 

and we obtain differential/difference properties of the orthogonal polynomials with respect 
to (2) in a unified way independently of the type of operator 3) , i.e., the results hold for 
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the continuous, discrete, and g-Hahn cases. We also find some relations between the 
semiclassical Sobolev orthogonal polynomials with respect to the inner product (2) and 
the semiclassical orthogonal polynomials with respect to u. To do this, we use recent 
results, obtained in [4,7]. 

In the last section of the paper we illustrate our results, applying them to some known fam- 
ilies of Sobolev orthogonal polynomials (Jacobi-Sobolev and A-Meixner-Sobolev ones) 
as well as to some new ones introduced in this paper for the first time such as g-Freud- 
Sobolev ones and another family related to a 1-singular semiclassical functional. 

In conclusion, we show that it is possible to unify the continuous, the discrete and the 
g-semiclassical Sobolev orthogonal polynomials by using a suitable notation, which we 
believe could be used in the future for further research. 



2 Basic definitions and notation. Semiclassical orthogonal polynomials 

As we said in the introduction, the Sobolev polynomials, the A-Sobolev polynomials 
and the g-Sobolev polynomials have been considered in the literature with different ap- 
proaches. Our main idea is to establish the main algebraic properties of the polynomials 
which are orthogonal with respect to the inner product defined in (2) in a gen- 

eral way where the linear functional u satisfies a distributional equation with polynomial 
coefficients: 

^{(f)u) = iPu, (3) 

where ^ is the differential or the difference or the g-difference operator. Notice that a 
functional u satisfying (3) with (p and polynomials, degV' > 1, is called a semiclassical 
functional. 

Since, depending on the case, the polynomials are orthogonal with respect to different 
linear functionals, we need to introduce a rigorous notation to embed all of them. 

(1) The Lattice. For the general case we use the variable z. Notice that in the g-Hahn 
case we need to replace z by x{s) = Ci(g)g* + C2(g)g~* + 03(5), where s = 0, 1, 2, . . . , 
with q e C, |g| 7^ 0, 1, and Ci, C2 and C3 could depend on q. However, in this paper 
we only consider the g-Hahn case, and so we work with the g-linear lattice 

x{s) := ci(g)g' + 03(5), 

where s = 0, 1, 2, . . . , and g G C, |g| 7^ 0, 1. 

(2) The operators. For the general case we use five basic operators: its dual the 
identity J^, the shift operator S''^, and For a better reading of the paper we 
collect these operators acting on a polynomial r G P in the following table: 
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(^+r)(x) 


(^-r)(x) 


Continuous case 


dr{x) 
dx 


dr{x) 
dx 


r(x) 


r(x) 


Discrete case 


(Ar)(x) 


(Vr)(x) 


r(x + 1) 


r(x — 1) 


g-Hahn case 




(^i/gr)(x) 


r{qx) 


r{x/q) 



where the backward and forward difference operators are defined as 

(Ar)(x) : = - J^)(r)(x) = r{x + 1) - r{x), 
(Vr)(x) : = (J^ - f?~)(r)(x) = r(x) - r(x - 1), 

and the g-difference operator is (see, for example, [5, Eq. (2.3)]) 

' r{qx) — r(x) _^ ^ 

r G P, g G C, |g| 7^ 0,1. 



(^,r)(x) 



(g — l)x 
r'(0), X = 0, 



(3) The constants. Throughout the paper we use some constants: g, [n], and [n]*. In the 
continuous and discrete cases we set g = 1, and [n] = [n]* = n, and in the g-Hahn 
case we set g G C, |g| 7^ 0, 1, [n] = (g" - l)/(g - 1), and [n]* = (g"" - l)/(g~^ - 1). 

Let P be the hnear space of polynomials and let P' be its algebraic dual space. We denote 
by (u,p) the duality bracket for u G P' and p G P, and by (u)„ = (u, x"), with n > 0, the 
canonical moments of u. 

Definition 2.1 For u G P', vr G P, and c G C, let vru, (x — c)~^u, and be the linear 
functional defined by 

(7ru,p) := (u,7rp) , p G P, 
p{x) — p{c) 



x-c) u,p : 



u. 



X — c 

{^*u,p) ■= -q{u,^p) , pe¥, 
and thus (^u,p) = —q^^ (u, !^*p). 



, PGP, 



(4) 



Notice that for any u G P', 
and 



[x — c){{x — c) ^u) 



u. 



(x - c) ^((x - c)u) = u - (u)o5c, 
where Sc is the Dirac delta functional defined by {Sc,p) := p(c), p G P. 

Taking these definitions into account we get for any u G P' and any polynomial n 

^(ttu) =(^+7r)^u + (^7r)u, 
^*{nu) =(^l'~7r)^*u+ (^*7r)u. 



(5) 
(6) 
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Furthermore, for any two polynomials, vr and ^, the following relations are fulfilled: 

^(7rO = (^+7r)^e + (^vr)e, (7) 
^*(7re) = (cf "7r)^*e + (^*vr)e. (8) 

Definition 2.2 Given an integer 5, a linear functional u, and a family of polynomials 
(Pn), we say that the polynomial sequence is quasi-orthogonal of order 6 with respect 
to u when the following properties are fulfilled: 

• If \n — m\ > 6 + 1, then {u.PnPm) = 0. 

• If \n — m\ = 6, then {u,pnPm) 7^ 0. 

Notice that if u is a semiclassical functional and is orthogonal with respect to u, i.e. 
(5 = 0, then is called semiclassical orthogonal polynomial sequence. 

Now we introduce some results regarding semiclassical orthogonal polynomials which allow 
us to obtain several results concerning Sobolev orthogonal polynomials in Sections 3 and 
4. There is a lot of literature on semiclassical polynomials and an authoritative paper on 
this topic is [14]. 

First of all, we define the concept of an admissible pair which appears frequently linked 
with the concept of the semiclassical linear functional. 

Definition 2.3 We say that the pair of polynomials {(piip) is an admissible pair if one of 

the following conditions is satisfied: 

• deg ip 7^ deg — 1 . 

• deg-i/^ = deg0 — 1, with 

ap + q-^[n\*bt^Q, n > 0, 
where ap and bt are the leading coefficients of and (p, respectively. 

We denote by a := max{deg(j) — 2,degip — 1} the order of the linear functional u with 
respect to the admissible pair (0, ^Z'). The class of such functional is the minimum of the 
order from among all the admissible pairs (see [7,13]). 

Remark 2.1 Taking into account the two possible situations for the admissibility we use 
the following notation to unify them. We write 

(j){z) = 6,+2;2"+' + ■ ■ ■ , V(^) = a^+iz"^' + ■ ■ ■ , 

where: 

• If t = a + 2 and p < a + 1, then we set 60-+2 = bt and a^+i = = ■ ■ ■ = a^+i = 0. 

• If p = cr + 1 and t < a + 2, then we set a^^i = ap and 60-+2 = ^cr+i = ■ ■ ■ = bt+i = 0. 

Definition 2.4 We say that the distributional equation (3) has an rig -singularity if there 
exists a non-negative integer Hq such that Ocr+i + q~^[no]*b„+2 = 0. Otherwise we say that 
the distributional equation is regular. 
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The following result will be very useful for our purposes. In the next result and in Lemma 
2.1 we give some details in their proofs for each case, namely the continuous, discrete and 
g-Hahn cases, with the objective of making them more readable. It is also important to 
point out that all the results presented throughout the paper are valid for all the cases 
with no additional restrictions. 

Theorem 2.1 Let u be a linear functional satisfying the distributional equation 

^{(j)u) = ipu, 

where cf) and ip are polynomials; then u also fulfills the distributional equation 

^*(0u) = V^u, 

where (f) and tp are polynomials, with deg(f) < a + 2, and degtp < o" + 1, with a = 
max{deg0 — 2, degip — 1}. 

Proof: The continuous case is trivial since 3i = 3i* = d/dx. 

In the discrete case ^ = A and ^* = V, so 

^(0u) = ^ A(0U) = +(0u) - 0U = +(0u) = {il) + 0)u. 

Observe that, due to the property of the shift operator , the latter expression is also 
equivalent to 0u = ^'^{{ip + 0)u). Thus, setting (p := (p + ip and ip = ip, we get 

^*(0U) = V(0U) = (0 + 1p)u - "((0 + 1p)u) =1pU = ^U. 

The g-Hahn case is a little more complicated, due to the notation, but the proof is almost 
identical. 

In this case ^ = A/Aa;(s), and ^* = V/Va;(s); thus the distributional equation for ^ 
can be written as 

^+(0u) (Pu {Ax)ijn ^ ^ ^+{(pu) {(P+{Ax)tp)^ 



Ax Ax Ax Ax Ax 



Then, 



{f-{l+{Axm (^'(0+(Ax)^)) (£-(0+(Axh^)) 

— Vu = — u — u 

Vx Vx Vx 

^ {^-{(P+jAx)^)) ^ _ 0u ^ /^_^ 

Vx Vx \ Vx / 

Thus, using the above expression, (6), and setting (p := (p + {Ax)ip, we get 



V (0+(Ax)^A)u / y,x v(d+(AxU) 

Vx \ Vx / Vx 

Therefore, it is a straightforward computation to check that both (p and ip are polynomials, 
with deg (p < a + 2, and deg ip < a + 1 where a = max{deg — 2, deg -0 — 1}. □ 
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Remark 2.2 • Although the reciprocal of this result is also true, and it can be proved 
in a straightforward way, we have decided not to include it since it is not used in this 
paper. 

• In the sequel, given a pair {(f), ip) we use the notation (p := (p + ip and ip := ijj for the 
discrete case, and (p := (p + {Ax)'ip and ip := qip for the q-Hahn case. 

The inner product associated with the semiclassical functional u, 

(p,r) := (u,pr) , 

is Hankel, i.e., 

{zp,r) = {p,zr), (9) 

and then the corresponding monic semiclassical orthogonal polynomial sequence (p„) sat- 
isfies the recurrence relation 

ZPn{z) = Pn+liz) + BnPn{z) + CnPn-l{z), n>0, (10) 

with initial conditions po = = 0. Moreover the square of the norm of p„, d"^ = (u,p^), 

satisfies the relation 

dl = CX-i, n>l. (11) 

It is well known that there are several ways to characterize semiclassical orthogonal poly- 
nomials. In our framework, we use a characterization of these polynomials that is usually 
called the first structure relation. 

Theorem 2.2 [12] & [7, Prop. 3.2] Let (pn) be a sequence of monic orthogonal polyno- 
mials with respect to u, and (p a polynomial of degree t. Then, the following statements 
are equivalent: 

(i) There exists a non-negative integer a such that 

n+t 

(p{z)&Pn+l{z) = J2 ^n,uPu{z), U > a, 
v=n—(j 

Xn,n-a 7^ 0, U > a + 1. 

(ii) There exists a polynomial ip, with degip = p > 1, such that 

^{(pu) = ipu, 

where the pair {(p, %p) is an admissible pair. 

We also use a recent result that establishes a nice relation between the semiclassical 
orthogonal polynomials (p„) and the polynomials (^p„+i). 

Theorem 2.3 [4, Th. 3.2] & [21, Th. 3.5] Let {pn) be a sequence of monic orthogonal 
polynomials with respect to u, and (p a polynomial of degree t. Then, the following state- 
ments are equivalent: 
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(i) There exist three non-negative integers, a, p, and r, with p>l,r>a + t + l, and 
a = max{t — 2,p — 1}, such that 

n+a n+cr 
Yl ^n,uPu{z) = J2 '^n,uPl}\z), 
u=n—a u=n—t 

where p^^\z) := [n + l]-^{^Pn+i){z) , ^n,n+a = ^n,n+a = I, u > max{o-, t + 1}, 

r,r—a^r,r—t 7^ Oj 

(^(0u),p^) = 0, p+l<m<2a + t + l, (^(0u),pp) ^ 0, 

and if p = t — 1, then limq^i(u,Pp)~^(u, 7^ mbt, m G N*, and where ht is the 

leading coefficient of (p (admissibility condition). 
(ii) There exists a polynomial ip, with degip = p > 1, such that 

S'{(j)u) = ipu, 

where the pair (0, ip) is an admissible pair. 

In this paper we use some linear functionals related to the semiclassical functional u. Now, 
we introduce them and give some properties that they satisfy. 

Lemma 2.1 Letu be a linear functional satisfying (3); then the linear functional {(pu) 
satisfies the following distributional equation: 

^(0^+(0u)) = (^0 + g^+V^)(f+(0u), 

setting q = 1 in both cases, continuous and discrete. 

Proof: The result follows by applying (5). In fact, in the continuous case, since S''^ is the 
identity, we get 

In the discrete case, S'~^ is the shift operator, so applying (5) we get 

^(0(f+(0u)) = A(0(f+(0u)) = ^+0 A(^+(0u)) + A0^+(0u) 

= + A0)^f+(0u). 

In the g-Hahn case, ff"*" is also the shift operator. Thus applying (5) we get 

Taking into account that in this case ^S''^ = qS'~^^, since the lattice x{s) is g-linear, we 
get 

^(0^+(0u)) = (^0^+(0u) + g(f+0(f+^(0u)) = (^0 + g(?+V^)^+(0u). 



8 



□ 

Lemma 2.2 Let u be a linear functional satisfying (3), and let \ he a linear functional 
such that u = S'^{(f)\); then v satisfies an analogous distributional equation. 

Proof: It follows from making computations in (3) where it is necessary to take into 
account how the operators & and S'^ act, that is, = q<S^Sl . In fact, we obtain 

□ 

So, the previous lemmas provide motivation for defining a family of linear functionals 
which will be used throughout the paper. 

Definition 2.5 Given a semiclassical functional u satisfying (3), fork E 1^ we define the 
linear functional in a constructive way as 

Uo = u, 00 = 4>, 

where (pk is a multiple of (pk-i- 

Remark 2.3 Notice that, by Lemma 2.1 and Lemma 2.2, when u is a semiclassical func- 
tional then Ufc is a semiclassical functional for every integer k. On the other hand, using 
Lemma 3.3 in [13] for the continuous case and Lemma 2.2 in [7] for discrete and q-Hahn 
cases, the statement of in terms of u^-i is well defined since if u is any semiclassical 
functional satisfying the distributional equations 

^(0iu) = %jjiu and ^(02u) = ^()2U, 

then for := gcd(0i,02) there exists a polynomial ip with deg-^ > 1 such that ^(0u) = 
■^u. Thus, observe that (pk is not determined in a unique way. 

As a consequence of Lemmas 2.1 and 2.2 and some straightforward computations we can 
deduce: 

Theorem 2.4 Say we are given a linear functional u satisfying (3), that is, 

^{(j)u) = ipu, 

where ip and (p are polynomials in the conditions of Definition 2. 3. Then, for any integer 
k and any polynomial tt we get 

^(vr^fcUfe) = TfUfc, 

where if is a polynomial of degree, at most, degn + ak + 1, (Jk being the order of the 
linear functional Uk with respect to the pair {(pki'^'k)- Furthermore, if the pair {(pki'^Pk) is 
admissible, then deg vf = deg tt + (7^ + 1 ■ 

Finally, we obtain a useful result for our purposes. 
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Proposition 2.1 // (pn) is a sequence of semiclassical polynomials orthogonal with re- 
spect to the linear functional u of order a, then the sequence of polynomials {^Pn+i) is 
quasi- orthogonal of order a with respect to the semiclassical functional w\. 



Proof: Using (3), (4) and (8), we have 

= -(u,p„+i(gV'^'" + = 0, 

when m + cr < 72. □ 

Remark 2.4 As far as we know, the reciprocal of Proposition 2.1 is no longer true without 
the assumption that v is quasi- definite where v is the operator with respect to which 
(S'Pn+i) is quasi-orthogonal. 

For every integer k we denote by (p^^-^) the sequence of monic orthogonal polynomials 
with respect to the linear functional u^. Thus, using Theorem 2.4 and Proposition 2.1, 
we have that the 2ip\^^ are quasi-orthogonal of order ak with respect to u^+i. Therefore, 
we get for every integer k, 

^pi+i= E c^n,k,uP?^'\ keZ. (12) 

u=n—af^ 



Then, using (12) and Theorem 2.3, we deduce for k = —1 the following relation: 

n+(T_i n+(T_i 

E cM~'^i')= E <Mz), (13) 

u=n--a-i u=n—t—\—cr—\ 

where ^*„_,_^_^„i = ^nn+a-x — 1; ^-i '■— deg0_i, and (J„i is the order of the functional 
u_i associated with the admissible pair (</)_i, ^/^.i). 

Observe that in the classical case, that is, cr = 0, using Theorem 2.4 we have (T_i = 0. 



3 The semiclassical Sobolev orthogonal polynomials 



Let us consider the ^-Sobolev inner product defined on P x P by 

{p,r)s = {vL^pr) + \{\i,Sip3ir), (14) 

where u is a semiclassical functional of order a and A > 0. We denote by (p„) the monic 
orthogonal polynomial sequence with respect to u and by (Q^^^) the OPS associated with 
the (^-Sobolev) inner product (-, ■)s which we call the semiclassical Sobolev orthogonal 
polynomial sequence. 
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Proposition 3.1 Let (pn) be a monic semiclassical orthogonal polynomial sequence and 
let (QI^^) be the semiclassical Sobolev orthogonal polynomial sequence. The following re- 
lation holds: 

n+(T-i n— l+CT-i 

E C,.Py'Hz) = Qi'l^A^) + E On,uQi'\z), n>a., + H\ (15) 

v=n—a-\ u=n—a-\—H'' 

where H* := max{t_i, cr_i}, 

(lo) 

with := (^QI^^{z)Q^^{z)'^ ^, and the other coefficients 6n^u in (15) can be computed 
recursively. 

Proof: If we apply (13) and expand it we get 

E c,.pI'^(')= E <,uPu{') = qL'I-A^) + E ^n.gS'H^), 

u=n—a—i v=n—t-\—(7-\ j=0 



where 



\u=n—a—i I \u=n—t—i—a—i 



Thus, using (12) and the orthogonality property we get for < i < n ~ — H* 



?-t + (T_l 



\ \yU=n—a-i j I 

a/u, I e" Wf^)w\ = o, 



, v=n—2a- 



and therefore (15) follows. To obtain (16) we only need to take i = n — a_i — H* in the 
above expression. 

Finally, if we apply this process for n — cr_i — H* < j < n + a-i — 1 we get an upper 
triangular linear system such that we can compute the coefficients Onj in a recursive way 
and hence the result holds. □ 

Next, we particularize the above result for the classical case, that is, we name as the 
classical Sobolev orthogonal polynomial sequence the semiclassical Sobolev orthogonal 
polynomial sequence where the linear functional involved has order a = 0, i.e., it is a 
classical functional. 
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Corollary 3.1 Let (pn) be a classical orthogonal polynomial sequence and let (Ql^^^) be 
the classical Sobolev orthogonal polynomial sequence. The following relation holds: 

plr'H^) = QL'\z) + fn-i{X)Qi\{z) + en-2{X)Qi%{z), n > 2, (17) 

with Q\'^\z) = pi{z), i = 0,1, and 

en-2(A) = e„^, eo(A) = e2, (18) 

/n-l(A) = {Sndl_i + iniSn-l - /n-2 ( A) )c/^„2 } > /o(A) = ^1, 

where 

Pi'^H^) = Pn{z) + 5nPn-l{z) + enPn-2{z). (19) 

Furthermore, in the classical case we can give nonlinear recurrence relations for the square 
of the norms of the Sobolev polynomials and for the coefficients e„(A) appearing in (17). 

Corollary 3.2 For n > 2, we have 

dl = dl+ [X[n]^ + 6n{Sn - fn~i{X)))dl_, + en{in " e„_2 (A) - (A) (<5„_i - /„_2(A)))d2_2, 

(20) 

with cIq = rfg, and = df + Xd^, and 

~ X[n]^ + Cn + U~Sn - /n-i(A)) + e„(e„ - e„_2(A) -/„_i(A)(5„„i - /„_2(A)))C„"ii ' 

(21) 

with the initial conditions 

eo(A) = 6~2, ei(A) = (22) 

Proof: Using (17) and (19), we get 

={Qk'Hz), Qk'Kz))^ = {Qi'\z),pt^}{z))^ = dl + X (u, (^QW)(^) {^Pi~'^){z)) 
=dl + A (u, (^gi-'))(^) {[n]pn-i{z) - ~6n{^pn-i){z) - e„(^p„_2)(^ 
=dl + X[n]'dl, + 5„ (u, gW(z) pn-i{z)) + e„ (u, QW(z) Pn-2{z)) 
=dl + X[nfdl_^ + 5„(5„-/„-i(A))dLi + e„(e„-e„_2(A)-/„_i(A)(5„„i-/„„2(A)))rf2 



n-2i 



which proves (20). Finally, (21) follows from (20) using (11), and the initial conditions 
(22) can be deduced from (18). □ 

The inner product (14) does not satisfy the Hankel property {zp,r)s = {p,zr)s, and 
therefore the polynomial sequence (Ql^^) does not fulfill in general a three-term recurrence 
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relation. However, we find an operator ^ which is symmetric with respect to the inner 
product (14), that is, (^p, t)s = (p, ^'r)s, where p, r G P. Thus, we generahze some 
results in this direction obtained for the continuous and discrete cases (see the surveys 
[10], [11], [15], [16], [18], and the references therein, and also [20]). 

Proposition 3.2 Let ^ he the linear operator 

J := {<g-^)J + -(^*0 - V^)^* - A(^-0)^ (23) 

where is the identity operator. Then, 

((cf-0)p,r)s = (u,p^r), p, reP. (24) 

Proof: According to Theorem 2.1 the linear functional u satisfies the distributional equa- 
tion ^*({/)u) = ?/^u. From (14), we have 



((^-0)p, r)s = (u, (^-0)pr) + A(u, ^((^^</.)pj ^r). 
Now, using relation (7) with vr = ^*(r) we get 

^((^-0)p^V) = ^(((f-0)p) + ((^"0)p)^(^*r), 

and so, 

{{^~^)p,r)s = (u,p (((f'0)r - A((f-0)^(^*r))) + A(u, ^((^-0)p^*r)). 

On the other hand, if we apply the relation (6) to the distributional equation with vr = 

we get 

= (^"0)^*u+ (^*0)u, 

and therefore 

(^^0)p,r)^ = (u,p((^"0)r-A(^~0)^(^*r)))--^(u,p(^/^-^*0)^*r) = (u,p^r). 



Corollary 3.3 The following identity holds: 

{{i/j - ^*(j))p,r)s = {^*u,p ^r), p,reF. 

Proof: Using the same properties as in the proof of the previous result, we get 

{{ijj - ^*^)p, r)s = (u, (V' - 3>*4>)pr) + A (u, ^((^ - 3>*4>)p) 



□ 



{^*u,p^r). 
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□ 

In the next result we prove that the hnear functional ^ plays for the inner product (14) 
a role equivalent to the one played by the multiplication operator by z for standard inner 
products (see (9)). 

Theorem 3.1 The linear operator ^ defined in (23) satisfies 

k/Vi r)s = {P, ^r)s r e P. 



Proof: Applying Proposition 3.2 and Corollary 3.3 together with (7) and the fact that 
{^P,r)s = ((^^0)p,r)5 --{{^- ^*0)(^»,r)s - A((^-0)^(^», r)^ 

q 

= {u,p^r) - ^(^*u, {^*p) ^r) - A(u, (^(^») ^ r) 



{u,p^r) + A(u, {&p){^^r)) 



□ 



Remark 3.1 Taking into account the proof of the last results, if gcd{S'~(j), ^*(f) — tp) = 
d > 1, then we can consider {S~(j))/d and {^*(f) — ip)/d in the definition (23) of the 
operator ^ and these results also hold. 

To end this section, we give some algebraic differential/difference results. One of them al- 
lows us to give an expression for the polynomial (^^{^z)pj^z) in terms of a finite number 
of ^-Sobolev orthogonal polynomials. Another one gives the second-order 3l -equation 
(differential or difference or g-difference equation) that is satisfied by the polynomials 

Corollary 3.4 The following relations hold 



i'=n+deg (p 

(^-0)(2)p„(^) = ^ fin,uQi'\^)^ ri>H, 

i/=n—H 
n+H 

/Q^n\z)= E ^^n,uPu{z), n>deg0 

v=n—dcg (f> 
n+H 

v=n~H 



where H := maxjdeg?/' — 1, deg0}. 
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Proof: To prove the first relation, we write 



n+deg (p 
i/=0 

where fJ^n,ud.l = {{^~ (p)Pn, Qi^^) s = {'^iPn^Qi^'^) which, using Proposition 3.2, vanishes 
when u + H < n since deg ^tt = H + deg ir. 

To prove the second relation it is enough to take (24) into account to get 

n+H 

where 'dn,vd1 = {\i,Pu ^ Q^^^) = {{S'~(j))pu,Qi^^)s which vanishes when u + degcj) < n. 
Finally, using Theorem 3.1 we obtain 

n+H 

where Wn^^d-l = {^Ql^\Qi'''^)s = {Qn \ ^Q^v^)s which vanishes when u + H < n, and 
thus the third relation is proved. □ 



4 The examples 



We illustrate the results obtained in this paper with several examples covering continuous, 
discrete and g-Hahn cases. Other examples appearing in the literature are also included 
in this general approach and their properties can be deduced from the results introduced 
here. We also remark that the nonstandard inner product associated with the g-Freud 
linear functional and another one related to a 1-singular semiclassical functional are new 
and they have not been considered before. 



4-1 A continuous example: the Jacobi-Sobolev Polynomials 



The family of monic Jacobi orthogonal polynomials (-P"''') is at the top of the continuous 
classical polynomials in the Askey scheme and they can be written as (see, for example, 

[8]) 



X] 



2"(a + l)n 
{a + P + n + l)n 



iFi 



-n, a + /3 + n + 1 
a + 1 



a, (3 > 
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In fact, this family satisfies the following orthogonality property: 

_ 2^"+"+^+^r(n + ma + n + l)r(/3 + n + l)r(« + /3 + n + 1) ^ 



a + /3 + 2n + l)(r(a + /3 + 2n + 1))2 
where the linear functional u"'^ has the following integral representation: 

(u"''^,p) = j^^P{x){l-xY{l + xfdx, 

and satisfies the distributional equation: 

^((1 - x2)u"''^) = (/3 - a - x(a + /3 + 2))u"'^, 
where, as we pointed out in Section 2, ^ = 3^* = so t = 2, p = 1, a = 0, and H = 2. 
Then, we can consider the Sobolev inner product defined by 

{f,g)s= t f{x)g{x){l-xni + xfdx + x[' f'ix)g'{x){l-xni + xfdx, 



where a, f3 > —1, and A > 0. This nonstandard inner product has been considered in a lot 
of articles, and we refer the reader to the surveys mentioned several times during the paper 
for more details. We denote by (Q^'^) the sequence of monic polynomials orthogonal with 
respect to the inner product {f,g)s, which are called monic Jacobi-Sobolev orthogonal 
polynomials. 

The results obtained in this paper allow us to recover some relations between Jacobi and 
Jacobi-Sobolev orthogonal polynomials: 



pa 1,13 1^^^ _ Q<^^P(^nf.'^ _|_ On]n-lQn-l{^) + ^n,n-2Qn-2 (^) ) 



n+1 
v=n—2 

Moreover, according to (23) we can define the linear functional 

^a,/3 ^ _ ^2)^ ^x{a-(3 + x{a + - A(l - x^)^. 

ax dx'^ ' 

and Corollary 3.4 yields 

n+1 



v=n—2 



n+1 
v=n—2 

Observe that the minus signs appear due to the factor (1 — in ^ . 
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4-2 A discrete example: the A-Meixner-Sobolev polynomials 



Monic Meixner orthogonal polynomials can be written as (see, for example, 

/ 



(3) 

M„(a;;/3,c) = -^^^2Fi 
c — IP 



V 



-n, —X 

/3 



1 



(3>0, < c < 1. 



In fact, this family satisfies the following orthogonality property: 



where the linear functional u^^ has the following representation: 



(/3)nc"n! 



x=Q 



X) 



Tix + r 



and it satisfies the distributional equations 



^*{{x + /3)u*^) = [x[l - -j+ (3jM 

^(xu^O = (x(c-l) + /3c) u^^ 
where ^ = A and ^* = V, so t = 1, p = 1, a = 0, and H = I. 

Now let us consider the nonstandard inner product defined by 



(25) 



where /3>0,0<c<l, and A > 0. This inner product is known in the literature as the A- 
Sobolev inner product. We denote by {Q^^\x] (3 , c)) the sequence of monic polynomials 
orthogonal with respect to (25), which are called monic A-Meixner-Sobolev orthogonal 
polynomials. 

Like in the Jacobi case, the results obtained in this paper allow us to recover some relations 
between the families of orthogonal polynomials (M„(a;; /3, c)) and {Q[^\x] (3, c)): 

M„(x; /3 - 1, c) = QW(x; /3, c) + (A; f3, c)Q[\{x- /3, c), 
{x - l)M„(x; (3, c) = Qi'Ux; /3, c) + fitlnQL'K^l c) + ftt%-iQnfi{x; c). 
Moreover we define the operator as 

= - 1)^ + A (1 + /3c - x{c - 1)) V - A(a; - 1)AV, 
and applying the results of Section 3 we obtain 

^''Qi'\x; (3, c) = M„+i(a;; /3, c) + K'^^M^x; (3, c) + <,_iM„_i(a;; /3, c). 
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4-3 A q example: q -Freud type polynomials 

The family of monic g-Freud polynomials, (Pn), satisfies the relation [4] 

(^P„)(x(s)) = [n]P„„i(x(s)) + a„P„„3(x(s)), n > 0, 

where x(s) = g^, with < g < 1, ^ = ^g, P„i = 0, Pq = 1, and Pi(x) = x. 

So by Theorem 2.2 we get = 1, a = 2, and t = 0; hence (P„) is orthogonal with 
respect to the linear functional u'''^ of class 2 which fulfills the distributional equations 

^{u"^) = iju'^^ , degV = 3, (26) 

^*((1 + x{q - lH')u'^^) = qipu'^^. 

Furthermore, these polynomials are symmetric (see [4]) and satisfy the three-term recur- 
rence relation 

xPn = Pn+1 + CnPn-1) Ti > 1. 

In fact, a straightforward computation shows a„ = A'(g)g~"'c„c„_ic„_2, and the sequence 
(c„) satisfies the nonlinear recurrence relation 

q[n]cn-i + K{q)q~"-^^CnCn-iCn-2 = [n- l]cn + -ft'(g)g"""-^c„+ic„c„_i, n>2, 

with initial condition cq = 0, and for the sake of simplicity we choose the initial conditions 
ci and C2 in such a way that 

cl + C1C2 = 1, 

and we have limg^i- K{q) = 4. Therefore, iplx) = —K{q)q~^x^ and the linear functional 
u^^ has the following integral representation: 

^ > ((g-l)ir(g)g-V";g^)oo 



where 



f{x)d,{x) = (1 - g) E + (1 - 9) E 

'1 7^ — n ? — n 



fc=0 fc=0 

In fact, this family satisfies the following orthogonality property: 



u''^, PnPm) — 2CiC2 ■ ■ ■ Cn^nA 



Now, to illustrate our results we can introduce the nonstandard inner product defined by 



(/'^)^= /'/(")^(") ((,-l)K(,|,-3,4..,4, d,{x) 



1 00 



1 - - ((g-l)ir(g)g-3g4.;g4)^ 



(27) 



A / {%f){x){%f){x) 
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We denote by (Qf) the sequence of monic polynomials orthogonal with respect to the 
inner product (27), which we call the monic g-Freud-Sobolev orthogonal polynomials. 

The theory developed in the previous sections allows us to link the g-Freud-Sobolev poly- 
nomials with the g-Freud polynomials. Taking into account the distributional equation 
(26) and Proposition 3.2 we define 

J^'>^ = (1 - (g - l)ir(g)g-V)^ + ^ir(g)g-V^i/, - A(l - {q - l)K{q)q-'x^)^,^i/g. 
In this case H = deg = 4 and the results in Section 3 can be rewritten as 

n+4 

(1 - (g - l)i^(g)g-V)P„(x) = ^ <uQf{x), 

u=n~4 

n+4 

u=n—4 
n+4 

u=n—4 

4-4 A 1-singular semiclassical polynomials of class 1 

The family considered for this example was studied by J.C Medem in [17]. Such a family 
of monic polynomials, (S'„), which we call Medem polynomials, is orthogonal with respect 
to the linear functional w which satisfies the distributional equation 

^(x^w) = (-x^ + 4)w, (28) 

where 2i = 2i* = t = 3, p = 2, and so a = 1, with initial condition (w)i = (w, x) = 0. 

ax 

Moreover a straightforward computation shows that a^+i + ba+2 = 0, i.e., w is 1-singular. 
Indeed, w is the symmetrized linear functional associated with the linear functional b'-"^-' 
(see [3] Chapter 1, Sections 8 and 9), i.e., (w)2n+i = and (w)2n = (b''~2-')„, for any 
n > 0. 

Notice that the linear functional b^"-* has the following integral representation: 

(b("\p\ = — / P{z)z''e-Uz, a > -2, 
\ ' 27rz JT 

and thus, after a straightforward calculation, we get that the linear functional w has the 
following integral representation: 

(w,P) = -i^ / P{z)z-^e-^dz, 
2m JT 

where T := {z E C : \z\ = 1} is the unit circle. 
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This family can be written in terms of the Bessel polynomials as follows [17]: 



S2n{x) =j-^-j-bL ^ V) = r2„x^-"e^(^*)" (x^'^-'e-i 



[n - 2Jn 



on 

S2n+i{x) =-, r-a:i?i"^^(x2) = r2.+ia:3-"e^(^*)" (x^^-^'e 

in — - L_ V 



1 



where is the Bessel polynomial of degree n with parameter a, and r„ 7^ are the 
corresponding normalization coefficients for n > 0. 

We can introduce the nonstandard inner product defined by 

{f.9)s = 7^ / /(^)^^-'e-^c?^ + A / J\zy{z)z-'e-l^dz. (29) 

We denote by ) the sequence of monic polynomials orthogonal with respect to the 
inner product (29), which we call the monic Medem-Sobolev orthogonal polynomials. 

Taking into account the distributional equation (28) and Proposition 3.2 we define 

= + A(4x2 - 4)-^ - Xx'^. 

ax dx^ 

In this case H = deg = 3 and the results in Section 3 can be rewritten as 



n+2 

x^Sn{x) = Q^_^_^{x) + ^ fij^j^Qj^{x), 

i^=n— 3 
n+2 

u=n—3 
n+2 

^'Q'n{x)=Q'n+3{^)+ E <uQu{x). 

u=n—3 
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Appendix. Proof of Theorem 2.4 

We are going to prove the result for /c > since the case k < is totally analogous to 
this one. When k = we know that (j) and ip are polynomials of degree, at most a + 2 and 
(J + 1, respectively, where a = ctq is the order of u with respect to the pair (0, ip). Then, 
using (5) we get for any monic polynomial vr, 

^(vr^uo) = (^7r)0uo + (cf +7r)^(0Uo) = ((^7r)0 + (^+7r)^)uo, 

where vr := (^7r)0 + {S^~^7i)ip is a polynomial of degree, at most, degTr + a + 1. 
Moreover, if the pair (0, ip) is admissible, then the leading coefficient is 

q"\aa+i + q~^['m]*ba+2) 7^ 0, with m = degvr. 

If we assume the result holds for every < k < K, let us prove it for k = K. Thus, the 
linear functionals satisfy the distributional equation ^(0A,.Ufc) = ipk^ik for Q < k < K , 
and then we can obtain 

= (^0/<_i)ux + +^k_i)^+(0a'-iUa'-i) 

Therefore taking into account Definition 2.5 and Remark 2.3, there exist polynomials 
(t)K,'4^K satisfying 

^((/.^ux) = V'i^ux, (30) 

and a monic polynomial such that (j)K-i = ^k4>k, and now, taking (30) into account 
the polynomials ipx and ipK-i fulfill the relation 

^{(i)K-l) + g(^+^K-l) = ^{^k)<Pk + ^^{^K)i^K, (31) 

and we can check that i/^k is a polynomial of degree, at most, ctk + 1- Furthermore, using 
(31) and working in the same way as in the case k = 0, we get 

^(tt^xUk) = ((^7r)0i^ + {^+n)i>K)uK. 

Clearly, vr := (^7r)0x + {<^~^'^)i^K is a polynomial of degree, at most, m + + 1- 
In fact, if the pair {(f) k, iPk) is admissible the leading coefficient of vf is 

g™(a<,,+i + g-i[m]*6.^+2)7^0, (32) 
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where the coefficients a^j^+i and 60-^^+2 can be obtained recursively 

(33) 

Notice that in the continuous and discrete cases q = 1 and [m]* = m. Then, using the 
admissibihty condition of Definition 2.3, the above expression, and (32) in a recursive way, 
we get that the leading coefficient is equal to 

a^K+i + = 0(7+1 + + ^((T^ + 2)j 6^+2 7^ 0, A; > 0. 

Therefore in the g-Hahn case there are infinitely many values of q for which the expression 
(32) is different from zero for any nonnegative integers m and k. Hence the result follows 
for k>0. 

Observe that by (33) if there exists an integer k such that the pair {(pk,i^k) is admissible, 
then the pair (0^, ipi) is admissible for every £ > k. 
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